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ABSTRACT 

Context. The [Cii] 158 pm fine structure line is one of the dominant cooling lines in star-forming active regions. Together with 
models of photon-dominated regions, the data is used to constrain the physical properties of the emitting regions, such as the density 
and the radiation field strength. According to the modeling, the [Cii] 158pm line integrated intensity compared to the CO emission 
is expected to be stronger in lower metallicity environments owing to lower dust shielding of the UV radiation, a trend that is also 
shown by spectral-unresolved observations. In the commonly assumed clumpy UV-penetrated cloud scenario, the models predict a 
[C ii] line profile similar to that of CO. However, recent spectral-resolved observations by Herschel/HIFI and SOFIA/GREAT (as well 
as the observations presented here) show that the velocity resolved line profile of the [C ii] emission is often very different from that 
of CO lines, indicating a more complex origin of the line emission including the dynamics of the source region. 

Aims. The Large Magellanic Cloud (LMC) provides an excellent opportunity to study in great detail the physics of the interstellar 
medium (ISM) in a low-metallicity environment by spatially resolving individual star-forming regions. The aim of our study is to 
investigate the physical properties of the star-forming ISM in the LMC by separating the origin of the emission lines spatially and 
spectrally. In this paper, we focus on the spectral characteristics and the origin of the emission lines, and the phases of carbon-bearing 
species in the N159 star-forming region in the LMC. 

Methods. We mapped a 4'x(3'^') region in N159 in [C ii] 158 pm and [N ii] 205 pm with the GREAT instrument on board SOEIA. 
We also observed CO(3-2), (4-3), (6-5), *^CO(3-2), and [Ci] ^Pi-^Pq and ^P 2 -^Pi with APEX. All spectra are velocity resolved. 
Results. The emission of all transitions observed shows a large variation in the line profiles across the map and in particular between 
the different species. At most positions the [Cii] emission line profile is substantially wider than that of CO and [Ci]. We estimated 
the fraction of the [C ii] integrated line emission that cannot be fitted by the CO line profile to be 20% around the CO cores, and up to 
50% at the area between the cores, indicating a gas component that has a much larger velocity dispersion than the ones probed by the 
CO and [C i] emission. We derived the relative contribution from C^, C, and CO to the column density in each velocity bin. The result 
clearly shows that the contribution from C^ dominates the velocity range far from the the velocities traced by the dense molecular 
gas. Spatially, the region located between the CO cores of N159 W and E has a higher fraction of C^ over the whole velocity range. 
We estimate the contribution of the ionized gas to the [C ii] emission using the ratio to the [N ii] emission, and find that the ionized 
gas contributes < 19% to the [C ii] emission at its peak position, and < 15% over the whole observed region. Using the integrated line 
intensities, we present the spatial distribution of /[cii]/7fir. 

Conclusions. This study demonstrates that the [C ii] emission in the LMC N159 region shows significantly different velocity profiles 
from that of CO and [C i] emissions, emphasizing the importance of velocity resolved observations in order to distinguish different 
cloud components. 

Key words. ISM: lines and bands - ISM: kinematics and dynamics - Galaxies: Magellanic Clouds - ISM: individual objects: N159 


1. Introduction 


The [Cii] 158/im fine structure transition is one of the 
dominant cooling lines in photon-dominated regions (P DRs; 
iTielens & Hollenbachl 1 198.51 : ISternberg & Dalgarndn~995ll . and 
it is used to constrain their physical properties. Because of its 
strength, it also has a great potential for tracing star formation 
activ ity in external galaxie s and for estimating the star formation 
rate (Ide Looze et alj|201ll and references therein). On the other 
hand, changes in the metallicity affect the abundance of the ma¬ 
jor coolants including carbon, the electron densities in PDRs, the 
dust abundance that determines the UV penetration, the photo¬ 
electric heating rate, and the H 2 formation rate, resulting in a 


larger layer in the low -metallicity enyironment dBolatto et al.l 
1 19991 : iRollig et'aP l2006l) . This metallicity dependence is also 
shown obseryationally as a two-order-of-magnitude difference 
in the intensity rati o of [Cii] and CO among galaxies with dif- 
ferent metallicities ( Modiizukietayi99^ Poglitsch et alJl9^ 
llsrael et ^119961 : iMadden et alJ 19971 201 ll) . In order to under- 
stand the extragalactic [Cii] emission, which is the superim¬ 
posed emission from different physical enyironments, our inter¬ 
pretation needs to bridge detailed Galactic studies to obserya- 
tions of marginally or non-spatially resolyed external galaxies. 
The Large Magellanic Cloud (LMC) is one of the most important 
targets b ecause its different metallicity from our Galaxy (~ 1/3 
of solar; ICarrera et al.ll200^ enables us to inyestigate the influ- 
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Table 1. Observation parameters. 


Line 

Erequency [GHz] 

Instrument 

7f 


HPBW" ["] 


[K] 

‘'CO(3-2) 

330.5879653 

ELASH+ 

0.95 

0.69 

19.0 

20-50 

0.06 

CO(3-2) 

345.7959899 

ELASH+ 

0.95 

0.69 

18.2 

30-60 

0.04 

CO(4-3) 

461.0407682 

ELASH+ 

0.95 

0.61 

13.6 

30-120 

0.11 

[Ci] 3Pi-3po 

492.1606510 

ELASH+ 

0.95 

0.6 

12.8 

100-250 

0.11 

CO(6-5) 

691.4730763 

CHAMP+ LEA 

0.95 

0.42 

8.8 

200-700", 500-3500^ 

0.08", 0.04^ 

[Ci] VFi 

809.3419700 

CHAMP+ HEA 

0.95 

0.38 

7.7 

100-500", 500-3500^ 

0.51", 0.14^ 

[Nil] 

1461.1338000 

GREAT LI 

0.97 

0.67 

18.3 

(50-200)s 

(0.10)s 

[Cii] 

1900.5369000 

GREAT L2 

0.97 

0.67 

14.1 

9-24 

0.30 


“ Forward efficiency. ^ Main beam efficiency. Half power beam width. 

Typical (> 80% of the pixels) effective integration time for one pixel in a 20" resolution map. For data with CHAMP^, contri¬ 
butions from different pixels are counted as an independent integration time, resulting in a significant decrease of feff at the edge, 
that is covered by a small number of pixels only. 

® N159E. ^N159W. ® In 50" resolution map. ^ Median of the baseline noise. 


ence of the metallicity on the thermal balance and chemical com¬ 
position in the ISM, and its proximity (50 kpc as a median dis¬ 
tance of 729 references; NASA/IPAAC Extragalactic Database) 
allows spatially resolved studies of the large star-forming regions 
visible in the EMC. 


In most cases, integrated line intensities of [Cii] and CO 
emission are used as diagnostics in external galaxies. How¬ 
ever, velocity resolved mapping observations of the [C ii] emis¬ 
sion in Galactic regions have recently been made with the Het¬ 
erodyne Instrument for the Ear-Infrared (HIEI) on Herschel 
and are being made with the German REceiver for Astron¬ 
omy at Terahertz Erequencies (GREAT) on board the Strato¬ 
spheric Observatory for Infrared Astronomy (SOEIA), and they 
show complex line profiles of the [C ii] emission, typically sub¬ 
stantially wider than expected from the [C i] and low- or mid- 
J CO emission profiles. The velocity information is thus es- 
sential in order to identify the origin of the [C ii] emission 
dPerez-Beaupuits et al.ll201^ ICarlhoffl2013h . and the significant 
difference in the velocity prohle of the [C ii] emission from the 
low- or mid-J CO emission indicates that a substantial fraction of 
the [C It] emission is accelerated relative to t he quiescent mate¬ 
rial, e.g., it undergoes ablation or is blown off dPedes et al| 20Tc ; 


Mookeriea et al.H2012l: Okada et al.ll 2012HSchneider et alj|2012 : 


Simon et^l2012t iPilleri et al.ll2Q12h . These results emphasize 

the need of velocity resolved observations to distinguish differ¬ 
ent velocity components with their different physical conditions. 


Before the HIEI and GREAT observations only the pioneer¬ 
ing heterodyne observations from the Kuiper Airbor ne Obser¬ 
vatory (KAO) were available. iBoreiko & Betzl dl99lh observed 
17 locations in the EMC and showed that the [Cii] emission is 
approximately 50% wider than the CO(2-l) line. However, in¬ 
vestigating its spatial variation and the line profile in detail re¬ 
quires mapping observations and a higher signal-to-noise ratio 
(S/N) than was available at that time; HIEI provided the nec¬ 
essary higher S/N, but the mapping ability was still limited and 
allowed only a few strip maps to be observed across a few promi¬ 
nent regions in the EMC. 


Since CO lines are observable from the ground, more ve¬ 
locity resolved observations have been done for them. In the 
N159 star-forming region, which we focus on in this paper, phys¬ 
ical properties of the c louds have been studied using low-J CO 
mapping observations (I.Iohansson et alJfT998l : lMinamidani et al.l 


120081 ; iMizuno et aT]|2010l:lMinamidani et al.ll201 Ih or using sin¬ 
gle point observations of CO(4-3), CO(7 -6), and [Ci] ^Pi-^Pq 
and ^P 2 -^Pi emissions (iPineda e t al.ll200 8h. bu t detaile d line pro¬ 
files have not been discussed. iBolatto et al.l (l2000li present a 
large-scale map and velocity resolved spectra of [Ci] ^Pi-^Pq, 
but with a spatial resolution of only 4', which is insufficient to 
resolve the individual source components spatially. 

The aim of our GREAT core project on the EMC and Small 
Magellanic Cloud (SMC) is to study physical properties by dis¬ 
tinguishing the origin of the emission lines using velocity pro¬ 
files in active regions in the EMC and SMC and to investigate 
the effect of the low-metallicity environment on the properties 
and evolution of the interstellar medium (ISM). In this paper, we 
present observations in N159, focusing on the velocity profiles 
and the origin of the emission lines, and present the column den¬ 
sity ratio between C^, C, and CO using simple assumptions. We 
do not deal with any elaborate PDR modeling or with the metal¬ 
licity effects. A full PDR model analysis will be presented in a 
follow-up paper (Okada et al. in prep.). N1 59 is located south 
of 30 Dor as a par t of the molecular ridge (ICohen et al.lll988t 
iMizuno et al. 20W). Three CO cores are identified (N159 W, E, 
and S: l.lohansson et alJl99^ . N159 S shows no evidence of star¬ 
forming activity, while W and E possess O-type st ars and Hit 
regio ns, where massive stars formed a few Myr ago (IChen et al.l 
1201 oh . The results of 30 Dor and N66 in the SMC are presented 
in a separate paper (Requena et al. in prep.). 


2. Observations and data reduction 

2.1. [Cii] 158pm and [Nii] 205pm observations with 
SOFIA/GREAT 

We mapped 4'x(3'^') in the N159 star-forming region in the 
EMC, observing the [Cii] emission at 1900.5369 GHz (158/im) 
using the Germ an REceiver for Astron omy at Terahertz Erequen¬ 
cies (GREATPt iHevminck et al.l2012h on boa rd the Stratospheric 
Observatory for Infrared Astronomy (SOEIA; [Young et al.l2012h 
in July 2013, as part of the open time and guaranteed time in 
cycle 1 observations. Observations were performed during four 

' GREAT is a development by the MPI fur Radioastronomie and the 
KOSMA / Universitat zu Kdln, in cooperation with the MPI fur Son- 
nensystemforschung and the DLR Institut fur Planetenforschung 
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flights of SOFIA’S first southern deployment to New Zealand and 
the L2 channel was tuned to [Cii]. During the first two flights, 
we observed in parallel [N it] 205 jum with the LI channel. The 
observations were made in total-power on-the-fly (OTF) mode, 
with 0.5 to 2 second integration time per dump, at 6" step size. 
The total flight time for N159 observations was about 4 hours 
and the total integration time on source was about 1 hour. The 
OFF position for the observations was located southwest of the 
observed region (05^'38"53^ 3 -69°46'29" (J200 0)), which is 
the same OFF position used in IPineda et al.l (l2008h . 

Pointing on SOFIA is done via a combination of fast iner¬ 
tial gyro stabilization of the telescope and regular corrections on 
the optical guide cameras. Experience shows that pointing and 
tracking is better than 2 arcsecs. The alignment between the tele¬ 
scope optical axis, i.e., the guide cameras, and the Foresight of 
the GREAT instrument is checked and verified at the beginning 
of each observing campaign after installation of the GREAT in¬ 
strument by observing a suitable planet, and was 1.5 arcsecs. 

The data from the extended bandwidth Eourier transform 
spectrometer (XEETS; 2.5 GHz bandwidth with 44 kHz res- 
olution) were ca librated by the standard GREAT pipeline 
(iGuan et alj|2012h . Details of the observations are listed in Ta¬ 
ble □ A polynomial of second order was used to fit the base¬ 
line. We checked the results obtained with polynomials of first 
to fourth order as a baseline, and found that a first-order poly¬ 
nomial results in higher baseline noise, fourth-order polynomi¬ 
als introduce artifacts, and second- and third-order polynomials 
make no significant difference. After baseline subtraction, the 
data were spectrally resampled to 1 km s ' resolution and spa¬ 
tially resampled to 20" resolution using the baseline noise as the 
weight. In the analysis of the very weak [N ii] emission, we re¬ 
sampled the data to 2 km s ' and 50" resolution to obtain a better 
S/N. Data reduction up to the sp ectral resampling was done by 
the CLASS software (|Petvl2005l) , and the spatial resampling and 
further analysis was done with IDL. 

2.2. CO and [Cl] observations with APEX 

In addition to the GREAT [C ii] and [N ii] mapping, we carried 
out complementary o bservations with the Atacama Pathfinder 
Experiment (APEX0; iGusten et al.l l2006h of CO(3-2), CO(4- 
3), '^CO(3-2), and [Cil ^P|-^Pn usin g the MPIfR heterodyne 
receiver PLASH+ (iKlein et alj |20l4 , and CO(6-5) and [C i] 
^P 2 -^Pi usi ng the Carbon Heter odyne Array of the MPIfR 
(CHAMP^: lKasemann et al1l2006h in June, July, September, and 
November 2012. The backends of PLASH+ are FETSs with the 
spectral resolution of 38 kHz and 76 kHz for the 345 GHz and 
460 GHz band, respectively. CHAMP^ has two hexagonal ar¬ 
rays each with 7 pixels, simultaneously operated for the 620- 
720 GHz and the 780-950 GHz frequency range, and the back¬ 
ends are XFFTSs (1.5 GHz bandwidth with 212 kHz resolution). 
All observations were made in the OTF mode, and the OFF posi¬ 
tion is the same as in the GREAT observations. Observational pa¬ 
rameters are listed in Table[T] Pointing was established on nearby 
R-Dor with CO lines, repeated every hour, and the pointing un- 
cer tainty is less than 2.5". 

iJohansson et al.l (Il998h mapped 1' x 1' around the brightest 
core of N159 W in CO(3-2) with the Swedish-ESO Submillime¬ 
tre Telescope (SEST), but it does not cover the w hole region of 
our [C ii] observations. iMinamidani et al.l (l2008h show a larger 


^ APEX is a collaboration between the Max-Planck-Institut fur Ra- 
dioastronomie, the European Southern Observatory, and the Onsala 
Space Observatory. 


CO(3-2) map including another prominent core of N159 E ob¬ 
tained with the Atacama Submillimeter Telescope Experiment 
(ASTE), but these are position switched observations with a grid 
of 10" to 40" and are not fully sampled. Therefore, we included 
CO(3-2) in our APEX observations in order to have sufficient 
complementary data. 

A linear baseline was fitted in tfie ranges 150 - 225 kms^' 
and 255 - 350 km s ' for ELASH+, and 100 - 220 km s ' and 
255 - 400 km s ' for CHAMP+. After discarding spectra with 
excess noise, the data were spectrally resampled to 1 km s ' and 
spatially to 20", identical to the resolution of the [C ii] data. As 
with the GREAT data, reduction up to the spectral resampling 
was performed using the CLASS software, and the spatial re¬ 
sampling and further analysis was done with IDL. 

3. Results and discussion 

3.1. Spatiai and veiocity structures 

Eigure [1] shows the integrated intensity maps of CO(3-2), (4- 
3), (6-5), '3C0(3-2), [Ci] ^Pi-^Pq and 3P2-^Pi, and [Cii], over- 
layed with the contours of IRAC 8 /urn from the Spitzer data 
archive. We estimate the uncertainty of the integrated intensity 
(cT/ [K kms^']) based on the baseline noise of the spectrum at 
each position of the resampled map (crjms, TablellJ, and the ve¬ 
locity range used for the integrated intensity (220-250 kms^'). 
In Pig. [T] pixels that have an integrated intensity smaller than 
3cr/ were masked. 

As shown in the following, the velocity structure strongly 
varies across the observed area, and from fine to line. Tfierefore, 
tfie integrated intensity map has to be combined with spectral 
information to interpret the spatial structure. Spectra at selected 
positions (marked in the CO(3-2) and [Cii] maps of Fig. [TJ are 
shown in Fig.|2l Fig. [3 shows the channel maps of the [Cii] and 
CO(3-2) emission; and the position-velocity (p-v) diagram of 
CO(3-2) and [C ii] along two selected cuts are shown in Fig. 0] 

In a first step of our analysis we con firmed tfiat our CO(3-2) 
ma p reproduces brigfi t cores shown by IJohansson et al.l (1 19981) 
and IMinamidani et ^ (l2008l) . AtN159 W (position 1 in Figs.[T] 
and|2] pi in the following), CO(3-2) and all other lines show a 
wider (~ 8 km s ' in FWHM) velocity profile compared to sur¬ 
rounding regions. Regions around pi (p2, p3, and p4) show nar¬ 
rower profiles with different central velocities, which may blend 
at pi. The integrated intensity peak of CO(3-2) is shifted south¬ 
west compared to the peak in IRAC 8 pm, and the IRAC 8 pm 
has a cometary shape facing northeast. The spatial distribution 
indicates a stratification from tfie nortfieast as an ionizing front 
to the molecular cloud toward the southwest. On the other hand, 
CO(3-2) in the northeast of N159 W does not completely van¬ 
ish, but shows a weak peak (p3); a very weak east-west bridge 
along IRAC 8pm is still seen in CO(3-2) (including p8). The cut 
(a) in the p-v diagram of CO(3-2) (Fig. Ell goes through p2, p3, 
p8, and p9, and it shows that the weak bridge connects clumps 
at p3 and p9 in the velocity as well, indicating that these clumps 
are not independent but rather connected with each other. The 
core south of N159 W (p4) is not strong in CO(3-2). From this 
core, a long ridge continues toward the south up to the edge 
of the map with a large velocity gradient (cut (b) in Fig. 01 and 
pl5). On the western side, there is another north-south ridge con¬ 
nected to N159 W. It contains a strong core in both CO(3-2) and 
IRAC 8 pm at the southern edge (pi4, the H 2 O maser source). 
At N159 E, the CO(3-2) peak (plO) shifts to the east compared 
to IRAC 8 pm, and the extension to the southeast is also the east 
side of the IRAC 8 pm ridge. The structure to the southwest from 
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Fig. 1. Integrated (220 to 250 kms^*) intensity maps (colors, 20" resolution) overlayed with the contours of IRAC 8yum. The last (eighth) panel 
shows a three-color image composed of CO(3-2) (blue), [Ci] ^Pi-^Po (green), and [Cii] (red). The red lines outline the observed area. Blue 
asterisks in the CO(3-2) and [Cii] maps mark the positions whose spectra are shown in Fig.|2] Open blue squares and re d triangles in the CO(4-3) 
map show the position of stars with a spectral type of 07 or earlier a nd 07 to Bl, respectively , from iFarina et al.l ll2009h . the filled blue star marks 
the pos ition of LMC X-1, th e blue asterisks are compact Hit regions dlndebetouw et al.ll2004l) . and the red diamond marks the position of an H 2 O 
maser dLazendic et alJl2003l . 


N159 E (p9) fills the gap in IRAC 8/im, and is connected to the 
bridge coming from N159 W. 

The spatial structure of prominent cores in the integrated in¬ 
tensity maps of CO(3-2), CO(4-3), CO(6-5), and '^€0(3-2) is 
similar, but the relative strength between cores looks different. 
For example, the relative strength of p4 to pi is larger in CO(6- 
5) than in CO(3-2), indicating a difference in the physical con¬ 


ditions in the different cores. Moreover, the spectral profile of 
CO(6-5) shows a clear difference with CO(3-2) at some positions 
because the different velocity components trace different clumps 
with different physical properties. For example at p7 (Fig. IS, the 
CO(4-3) and (3-2) emissions have two significant velocity com¬ 
ponents at ~ 235 km s ' and ~ 240 km s ', while the CO(6-5) 
emission at ~ 235 km s“' just shows a wing, and its ~ 240 km s“' 
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Fig. 1. (continued) 


component is very strong in terms of the ratio to CO(4-3) or (3- 
2). At p6, both the CO(6-5) and [Ci] 3P2-^P 1 emission have a 
peak at ~ 237 kms“^ and especially the CO(6-5) emission is 
strong, while CO(4-3) and (3-2) have a broad velocity prohle 
peaking at > 240 kms^', and [Ci] ^Pi-^Pq is even not detected. 
Position 6 is close to LMC X-1, which is t he X-ray binary inc lud- 
ing a black hole with a mass of ~ 1 IMp (Orosz et al.ll200^ and 
an 07 III companion (ICowlev et al.lll995h . The ~~237 km s ' 
component may originate in a clump that is affected by strong 
X-ray emission from LMC X-1. 

The [C i] ^Pi-^Pfl map and spectra show an overall similarity 
with CO(3-2), but the two north-south ridges south of N159 W 


are stronger in [Ci] ^Pi-^Pq. Position 15 is the peak in [Ci] ^Pi- 
^Po, whereas the peak there is not clearly seen in other emission 
lines. The peak of the western ridge shifts slightly to the north 
in [C i] ^Pi -^Pq. The [C i] ^P 2 -^Pi emission is detected in a wide 
area of N159 W and around the peak of [Ci] ^Pi-^Pq in N159 E. 
It is detected in only a small area around N159 E, which may be 
because of the shorter observing time and thus the lower sensi¬ 
tivity there (Table [T]i. 

The [C It] emission shows significantly different velocity pro- 
hles and integrated intensity maps from the other tracers. It typ¬ 
ically has a wider velocity profile, except for a few positions. 
At pi and pl4, all the lines have a wider and narrower veloc- 
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Fig. 2. Spectra at selected positions marked in the CO(3-2) panel in Fig.[T] 


ity profile than the other positions, respectively, and there is no 
significant difference in the width between [C ii] and other emis¬ 
sions, except that the peak velocity of the [C ii] emission at pi is 
shifted to a smaller velocity (Figs.|2]and|4|i. At p9-pl2, the [C ii] 
emission is not only wider than in the other tracers, but the peak 
velocity is also shifted by a few km s '. The integrated inten¬ 
sity of [Cii] spatially matches the IRAC 8/rm; in N159 W, their 
peaks match around pi while those of CO and [C i] emission are 
shifted southwest, and in N159 E their peaks are instead at pi 1 
in contrast to the CO and [C i] peak at plO. An exception is pl3, 
where the [C ii] emission indicates a clump but there is no corre¬ 
sponding IRAC 8/im emission. All C O and [Ci] emissi on lines 
are also very weak, and the Ha map bvIChen et al.l(l2010h a nd the 
843 MHz radio continuum map by iMills & Turtle! jl9^ show 
the presence of ionized gas position, so that it is likely that the 
[Cii] emission at pi3 also originates in the ionized gas compo¬ 
nent. Unfortunately, our [N ii] map, which would allow the [C ii] 
contribution from ionized gas to be estimated, does not cover 
this position. The ridge-like structure between N159 W and E 
(including p8) is significant in [C ii]. The p-v diagram of cut (a) 
in Eig.|4]shows that the velocity gradient and the central velocity 
are similar between CO(3-2) and [C ii], indicating that the emis¬ 
sion along the cut in the two tracers does not come from fully 
independent clouds but is related, although the [C ii] emission 
has a much wider line width. At p7, the velocity component at 
~ 240 km s ', which is significant in CO(6-5), is not detectable 
in the [C ii] emission. 


We compare our [C ii] observatio ns with prev i ous st udies and 
archival data from Herschel/PACS. Ilsrael et ^ (Il996h mapped 
the N159 region with EIEI/KAO; they quote a beam solid an¬ 
gle as a 68"disk, which corresponds to a Gaussian EWHM of 
64". We convolved our map to 64" resolution and derive a 
line integrated intensity of (6.0-6.3) x 10“"' erg cm“^ s“' sr“' 
at the peak in N159 E and W, whereas their result w as (3.5- 
4.0) X 10“"' erg cm“^ s ' sr“'. iBoreiko & (Il99lh detected 
the [Cii] emission at N 159 EIR(NE) with their far-infrared het¬ 
erodyne Schottky-mixer receiver on the KAO, showing a peak 
antenna temperature of 5.3 K and an integrated intensity of 
4.8 X 10“"' erg cm“^ s“' sr“' in the 43" EWHM beam. Convolv¬ 
ing our map to 43" and extracting the spectrum at their mea¬ 
sured position (which is 25" off from our peak of N159 E), we 
obtain a peak main beam temperature of 5.9 K and an integrated 
intensity of 5.9 x 10“"' erg cm“^ s“' sr“'. The value measured 
with GREAT is 10-20% higher, but considering that we are com¬ 
paring the antenna temperature with the main beam temperature 
for an extended source, the two observations are compatible. We 
also retrieved [C ii] observations with PACS from the Herschel 
data archive, and averaged the spectra in a 20" diameter circle 
based on the geometrical overlap of the PACS spaxels. The val¬ 
ues derived for N159 W (pi and p4) and E (plO and pll) are 
(6.1, 7.7, 7.4, and 8.2) x 10“"' erg cm“2 s“' sr“' for PACS and 
(7.2, 10.8, 8.4, and 10.3) x 10“"' erg cm“2 s“' sr“' for GREAT, 
showing that the PACS intensity is 0.7-0.9 of the GREAT in¬ 
tensity. We consider this difference to be within the calibration 
uncertainties of both instruments. 
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Fig. 3. 1 km s ' wide channel maps of the [C ii] emission overlayed with the contours of the CO(3-2) emission. The contours start from 1.5 K with 
a spacing of 1.5 K. The central velocity is written in each panel. 


3.2. Quantitative anaiysis of the veiocity profiies 

To investigate quantitatively the velocity profiles and their differ¬ 
ences in different emission lines, we fit the CO(3-2) spectra with 
two Gaussians, and apply the fitted center velocity and width 
in the fit to the other emission lines, fitting only their ampli¬ 
tudes. Across the map, there is no position that has more than 
two clearly separate velocity components. In order to handle po¬ 
sitions where a de-composition into two velocity components is 
not possible, we first fit the CO(3-2) spectra with two Gaussians 
with center, width, and amplitude as free parameters, and when 
the separation of the fitted center of the two Gaussians is less 
than 1 km s ', or the fitted amplitude of one of the Gaussians is 
less than the noise of the baseline, we re-fit the spectra with a 
single Gaussian. For the positions where the two velocity com¬ 
ponents are not clearly separated, the fitting may not be unique, 
but the sum of the integrated intensity of the two components is 
robust. Once the width(s) and center(s) of the Gaussian(s) from 
the CO(3-2) fit are determined, we use these values as fixed pa¬ 
rameters in the fitting of the other emission lines, thereby leav¬ 
ing their amplitude as the only free parameter. In these fits, we 
restricted the spectral range to be fitted to half of the FWHM 
of the CO 3-2 fit, as long as there is a data point in this range 
(otherwise we expand the fitting range to the FWHM) to make 
sure that the peak of the line is well fitted when the profile of 
the line is different from that of CO(3-2). This procedure allows 
us to identify how well the other emission lines match the ve¬ 
locity structure visible in CO(3-2) and to identify whether the 


integrated line emission from the other tracers shows evidence 
for additional velocity components not visible in CO(3-2). 

Figure |5] shows the relation between the integrated intensity 
and the sum of the fitted Gaussians for each emission line. For 
CO(3-2), both estimates match perfectly, which guarantees that 
the fitting with two Gaussians is appropriate. For CO(4-3) and 
CO(6-5), the integrated intensity is also very well reproduced by 
the sum of the fitted Gaussians, confirming that the three CO 
emission lines contain the same two velocity components and 
give no indication of additional ones. On the other hand, the 
plots for *^CO(3-2) and [Ci] ^Pi-^Pq show that the sum of the 
fitted Gaussian provides larger values by 5-10% than the inte¬ 
grated intensity where the emissions are strong. The [C i] ^P 2 - 
^Pi emission may have the same trend for the six strongest data 
points, but the S/N of the line is much smaller and the trend is 
not conclusive. Looking at individual spectra, we find that this 
is because the line widths of *^CO(3-2) and [Ci] ^Pi-^Po are 
slightly narrower at these positions than those of CO(3-2). When 
we fit the spectra with the width as a free parameter, but fixing 
the center and the amplitude to the values that we get in the pre¬ 
vious fit, the integrated intensity is well produced in both lines. 
This is likely an optical depth effect because the CO(3-2) and 
CO(4-3) are optically thick at the line center, the optical depth 
of the CO(6-5) emission at the strongest position is about unity, 
whereas the *^CO(3-2) and [C i] emissions are optically thin (see 
Sect. 13.31 1. In any case, the difference in the line width is small. 
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Fig. 4. (lower left) CO(3-2) integrated intensity map overlayed with the contours of IRAC 8 pm (same as Fig. [3, together with two lines marking 
the two cuts, along which the p-v diagrams are shown in the upper left and lower right panel, (upper left / lower right) The p-v diagrams of the 
[Cii] (color) and CO(3-2) (contour) in [K]. The contour spacing is 1 K. The position axis is horizontal in the upper left panel (cut a) and 
vertical in the lower right panel (cut h), so that the position aligns with the position of the lower left panel. 


since the dynamics makes the line width broader than the ther¬ 
mal width. 

The largest difference is for the case of [Cii]. The plot for 
[C It] in Fig. |5] shows a clear trend; the sum of the fitted CO- 
constrained Gaussians only contributed a fraction of the total 
flux calculated by the integration over the whole velocity range. 
The large scatter in the plot also indicates a large variation in the 
fraction of the [C ii] flux in the velocity components correlated 
with the CO and [C i] emission. For example, at the center of 
N159 W (pi, see Figs. [T] and |2ll, ~ 20% of the [Cii] integrated 
line emission cannot be contributed by CO-constrained Gaus¬ 
sians, whereas this ratio is ~ 40% at p7 and more than 50% at 
p9. The origin for this [C it] emission may be due to gas ablating 
from the dense cloud cores, or may be due to additional gas not 
immediately associated with the dense cores. 

We also examined the ratio of the integrated intensity of the 
different emission lines derived as the sum of the fitted Gaus¬ 
sians to that of CO(3-2). They show clear spatial variations, 
which is already visible from the example spectra (Fig. la. The 
[Cii]/CO(3-2) ratio varies most significantly: more than one or¬ 
der of magnitude across the map. The ratio is low where the 
CO(3-2) emission is strong (0.6-1 forpl, plO , and pl4), indicat¬ 
ing a low radiation fie ld and/or a high density (iRollig et al.l2006l: 
iKaufman et al.l [200^ . and it increases towards regions where 
diffuse gas is likely dominant (~ 11 at p7). A full analysis of 


these amplitude ratios using PDR modeling will be presented in 
the follow-up paper (Okada et al. in preparation). 

3.3. Column densities of C+, C, and CO 

In this section, we present and discuss the spectral and spatial 
variations of C^, C, and CO column densities and their ratios. 
We derive column densities under the simplifying assumption of 
local thermodynamic equilibrium (LTE) for each species, and a 
uniform excitation temperature (Tex) for different velocity bins 
at each position. 

Assuming that CO and '^CO have the same Tex and beam 
(area) filling factor ( 77 ), the peak temperature (Tb) ratio of CO(3- 
2) and '^CO(3-2) is described by [1 -exp(-Ti2)]/[l -exp(-Ti3)], 
where T12 and T13 are the peak optical depth of the CO(3-2) 
and '^CO(3-2) emission, respectively. We assume an isotope 
ratio *^C/*^C of 4 9 as derived for the LMC/N113 region by 
IWang et ^(l 2009 h and use it as T12/T13. We derive T13 by solving 
Eq. ([T]i numerically pixel by pixel, 

rB(CO( 3 - 2 )) ^ l-exp(-49xTi3) . 

Tb) 1300(3-2)) l-exp(-Ti3) ’ 

ti 3 varies across the map from 0.06 to 0.5, which corresponds to 
ti 2 = 3-35, i.e., CO(3-2) emission is optically thick everywhere. 
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The optical depth of an emission line Ty is related to the col¬ 
umn density of the lower state Ni by 


- ^ Suj., . 

8 ^ Vq gi 


1 I 


( 3 ) 


where subscripts u and I indicate the upper and lower state, g 
is the statistical weight, Aui is the Einstein coe fficient, and (f>y is 
the normalized line profile (IWilson et al.ll2013l) . Under the LTE 
assumption 


A 

Mot 



( 4 ) 


where Mot is the total column density, £; is the energy of the 
lower state, and Z is the partition function 


^ 77 


( 5 ) 


Therefore, the column density in each velocity bin (dNy) can be 
derived as 


dNv — Tv X 


87 rv 2 z 1 ( E, \ 

guAui^^'^\kT^y] 


1 - exp 


hvQ \ ' 

kT,J 


dv, ( 6 ) 


Ty is derived from the brightness temperature at each velocity bin 
(7’b(v)) as 


Ty = - In 


1 - 


kTB(v) 

rjhv 




( 7 ) 


Fig. 5. The integrated intensity as a sum of the fitted Gaussians vs. the 
integration over 220-250 km s“*. The blue line corresponds to both val¬ 
ues being equal. 


In general, Tex can be expressed as 


hv 

T - _ 

^ ex , 

k 




-1 


( 2 ) 


Here we use T 13 and 7 'b(^^CO( 3-2)) to derive Tgx of CO. The 
resulting T^x for the case of rj - 0.5 is shown in Fig. [ 6 ^. In the 
figure, the position s of early-ty pe stars (earlier th an B4) from 
iF^ina et’afl (l2009t). YSOs fromjChen et al.l (l2010t) , and OB star 
candi dates from lNakaiima et al.l (12005 ) (taken from IChen et al.l 
I2OIOI1 are oyerplotted. Most of the cores are associated with one 
or more of these, which is consistent with a higher Tex at those 
positions. The calculation with 77 = 1 giyes a lower limit of Tex, 
and the Tex at N159 W is 21 K, 36 K, and 15() K for 77 = 1, 
0.5, and 0.1, respectiyely. iGalametz et alJ (1201 3h deriye a dust 
temperature from the far-infrared continuum emission with 36" 
resolution, which is 30-35 K in our obseryed region. When we 
estimate Tex at N159 W and E in the same spatial resolution, 
77 = 0.3-0. 4 giyes the same tempera ture for Tex as their dust 
temperature. iMinamidani et al.l (1201 ill deriye the kinetic temper¬ 
ature from a LVG analysis with 45" resolution to be 40-82 K 
and 65-151 K in N159 W and E, respectiyely. If we assume that 
Tex equals their kinetic temperature, it requires a yalue for 77 of 
< 0.3. The deriyed range of filling factors is consistent with the 
reasonable scenario, that clumps with an enhanced density by 
more than one or two orders of magnitude fill ~ 60% or ~ 35% 
of the projected area, respectiyely. In the following, we use the 
results for 77 = 0.5, together with 77 = 0.1 and 0.2 for some cases. 


In the case of the CO 7 -(- 1 —> 7 transition, Ei - hBJ(J - 1 - 1), 
gu = 2(7 + 1) -H 1 , Z = kTy^JhB , and = 647r"^/(3/7c^) x Vg/ 7 ^ x 
(7 + l)/(27 -H 3) (I Wilson et al.ll2013h . therefore 


dNy{CO) — Tv X 


3ckn 


1 


^TT^voBfE- 7+1 


exp 


hBJiJ + 1) 

kT„ 


1 I 


dv. 


( 8 ) 


where /i is the dipole moment and B is the rotational con¬ 
stant. We used rB(v) of '^CO(3-2) and deriyed 7My('^CO) and 
conyerted it to dNv(CO) by multiplying '^C/'^C= 49. The to¬ 
tal column density N(CO) was deriyed by integrating 7My(CO) 
oyer the frequency range corresponding on the yelocity of 220 - 
250 km s ' (Fig. [8]|. The spatial distribution of M(CO) follows 
the integrated intensity of CO(3-2) and (Fig.| 6 ^). 

In Fig.|2]we compare the obseryed peak temperatures of 
CO(4-3) and CO(6-5) with the expected ones from the peak (in 
terms of yelocity) of dNy{CO) assuming that all transitions haye 
the same Tex and 77 . In the case of 77 = 0.5, the obseryed Tb is 
well reproduced, indicating that the assumption of identical Tex 
and 77 for all transitions is justified. When 77 = 0.1, the obseryed 
Tb of CO(4-3) and CO(6-5) is systematically higher by ~ 20% 
and ~ 60%, respectiyely. The fact that the obseryed of CO( 6 - 
5) is weaker is likely due to CO(6-5) not being thermalized. The 
critical densities of CO(3-2) and CO(6-5) are (3^)xl0^ cm^^ 
and 3 x 10^ cm~^, respectiyely, for 10-100 K (calculat ed using 
lFlowed200lHGooryitchl 199^ . iMinamidani et ^(1201 ll) deriyed 
7 2(H?) of (4-5)xl0^ cm“^ for N159 W and E, and Pineda et~an 
(l2008h discuss the result of PDR modeling with an ayerage den¬ 
sity of the clump ensemble of about 10^ cm“^. These densi¬ 
ties are consistent with the assumption that the CO(6-5) emis¬ 
sion does not reach LTE, in contrast to CO(3-2). Calculations 
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Fig. 6. (a) CO excitation temp erature assumi ng a beam filling factor rj = 0.5 overlayed with contours of CO(3-2) integra ted intensity. The 
positions of YSOs identified bv IChen et al.|l|20 10|l are marke d with asteri sks, OB stars (earl ier than B4) from iFarina et al.l (l2009h are marked with 
triangles, and OB star candidates from lNakaiima et al.ld20Q^ . taken from IChen et in @13), are marked with circles, (b) Tex derived from the two 
[C i] lines overlayed with the contours of [C i] ^P| -^Po integrated intensity, (c) Pex derived by assuming t = \ for the [C ii] emission overlayed with 
the contours of [C ii] integrated intensity. 
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Fig. 7. Peak brightness temperature of CO(4-3) (left panels) and CO( 6 - 
5) (right panels) estimated from the CO column density derived by us¬ 
ing Tex under the assumption of 77 = 0.5 (upper panels) and 0.1 (lower 
panels) for all CO lines against observed Tg (20" resolution). Blue lines 
indicate the case where the estimated and observed Tg is the same. 


with the RADEX radiative transfer program show that the CO(3- 
2)/CO(6-5) intensity ratio at 30-100 K changes by more tha n 
a factor of 2 between 10^-10^ cm^^ (Ivan der Tak et al.ll2007h . 
which can explain the CO(6-5) discrepancy in Fig.|2]in this work. 
Alternatively, it is also possible that tj of CO(6-5) is smaller than 
that of CO(3-2) because it traces warmer gas which is expected 
to be in a thinner layer. 

At positions where both [C i] lines are detected, we derive 
the [C i] excitation temperature Tex and the peak optical depth r 
from their ratio and their peak temperature. The optical depth ra¬ 
tio of both transitions in [C i] can be expressed in terms of Tex un¬ 
der the LTE assumption (Eq.O. Combining with Eq.|2]for both 
transitions provides two independent equations in terms of Tex, 
77 , and Tio (the optical depth of ^Pi-^Pq) as unknown variables, 
and Tb as observed quantities. Since the S/N of [Ci]^P 2 -^Pi is 
low, and the line profiles of the two [C i] lines are the same with 
the exception of a few positions mentioned in Sect. 13.11 we did 


not use directly ^^(^Pi-^Po) and rB(^P 2 -^Pi) as observed vari¬ 
ables, but used ^^(^Pi-^Po) and /( 3 p 2 - 3 Pi)//(^Pi-^Po), where I 
is the integrated intensity. We calculated tiq and Tex numerically 
for individual rj. The resultant tiq is less than 1 at all pixels for 
all rj, and Tex is insensitive to rj. For the fully optically thin case, 
/(^P 2 -^Pi)//(^Pi-^Po) is proportional to T 21 /T 10 , therefore inde¬ 
pendent on rj. The median of Tex across the map is 70-80 K for 
77 > 0.2 and 90 K for 77 = 0.1. In Fig.| 6 ]r the spatial distribution 
of Tex for 77 = 0.5 is shown. In contrast to Tex derived from the 
CO emissions, Tex for [C i] is lower at the N159 W core than the 
surroundings. If we assume that Tex of CO and [C i] are simi¬ 
lar (within 10 K), the comparison at pi (N159 W), p4, and plO 
(N159 E) gives 77 of 0.3-0.4, < 0.2, and 0.1-0.2, respectively. 

To derive the total neutral atomic carbon column density 
from Eqs. (0 and (|7]i, we sum over the three fine structure levels 

Z = » + *,exp(-|L)+ 5 ,„p(-j|-) (9) 

Since the region where the [Ci] ^P 2 -^Pi emission is detected 
is small and the uncertainty of the derived temperature is rela¬ 
tively large because of the weak [Ci] ^P 2 -^Pi emission, we ex¬ 
amine both cases using the T^x derived from the [C i] ratio and a 
constant value of its median to estimate the column density. The 
derived column densities match within 10%, even less than 5% 
at most positions. In the following, we use the case with constant 
Tex (Fig.© because it covers a larger spatial area. 

Since there is only a single [C ii] line, we need to assume its 
optical depth in order to estimate Tex using Eq. I©, or assume a 
constant Tex for the column density derivation. In several Galac¬ 
tic PDRs, the optical depth of the [C ii] line is estimated to be 
(0.6-3.21 using the hyperfine components of the 11 I emis sion 
dOssenkopf et al.ll2013t iGraf et al.ll2012L IStacev et al.l|l991l) . In 

the N159 region, the strongest [*^Cii] hyperfine component is 
blended with the broad main isotope line, and the upper limit 
of the second strongest [*^Cii] hyperfine component does not 
give a useful constraint to the optical depth of [C ii]. Figure | 6 }; 
shows the estimate o f Tex for a constant optical depth t = 1 (see 
iKaufman et al.l 19^ for the peak velocity and 77 = 0.5. The spa¬ 
tial distribution follows the distribution of the [C ii] brightness 
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Table 2. Total column density of C, and CO, and Nu estimated from 
N{C^)+N{C)+N{CO). 
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because a stronger peak temperature is attributed to a higher tem¬ 
perature. When we assume r = 0.5, Tgx increases by 15-37%. 
For the optically thin case, t is coupled with rj. When we assume 
T - 3 and r » 1, T^x goes down typically by 10-16% and 11- 
18%, respectively. The latter case gives the lower limit of Tex. 

The column density of is derived using Eqs. ® and Qi, 
together with the first two terms in Eq. (|3 because is a two- 
level system. We examine both cases using the T^x derived by 
assuming t = 1 and a constant value (Eigs.[ 8 }; and d). The con¬ 
stant value is determined, depending on 77 , as the largest value 
of Tex estimated from the t » 1 case across the map, which is 
67 K and 199 K for j] = 0.5 and 0.1, respectively. In this way, the 
constant T^x is above the Tex lower limit at any position of the 
map. We then estimate t in case of this constant Tex. It is fully 
optically thick at the peak of p4, and regions around other 
peaks (N159 W, E, p7) have r ~ 1 at the peak velocity. Eor the 
constant T^x case, the contrast of the [C ii] emission is attributed 
to the difference in r; therefore, the derived column density ba¬ 
sically reproduces the [Cii] intensity distribution (Eig. |3l). Eor 
the constant t case, the stronger [C ii] emission is attributed to 
a higher Tgx (Eig. |^), but the column density is insensitive to 
Tex- Thus, the large variation of the [Cii] velocity profile causes 
a different Ty profile from position to position and results in a 
different spatial distribution of A^(C^) from the T^x and the in¬ 
tegrated intensity distribution. The higher A^(C^) in the eastern 
part of the map (Eig. | 8 };) is due to a wider line profile there be¬ 
cause the constant t is at the peak of the line and a wider line 
makes the integrated column density larger. Qualitatively, posi¬ 
tions with strong [C ii] emission are the outer part of PDRs where 
the temperature is higher and the column density of is higher 
than toward the molecular core. Therefore, the assumptions of 
constant r and constant T^x are the two extreme cases, and the 
reality should be somewhere in between. In the following anal¬ 
ysis we use the constant T^x case because it is unlikely that the 
optical depth is the same between cores or the [C ii] peak posi¬ 
tions and regions in between where the [C ii] emission is much 
fainter. 

Assuming that N(C^)+N(C)+N(CO) is the total column den¬ 
sity of the carbon species in the gas phase, we estimate Ah using 


Fig. 8. The total column density of CO, C, and for the rj = 0.5 case 
overlayed with the contours of the integrated intensity of CO(3-2) in (a), 
[C i] -^Po in (b), and [C ii] in (c) and (d). N(C) is derived by assuming 
Tex= 74 K, and for A(C^), both cases of t = 1 (c) and 7’ex= 67 K (d) 
are shown (see text). 



230 235 240 245 

Velocity [km s '] 


Fig. 11. Spectra of the fraction of CO, C, and column density 
against the total column density, i.e., dNy(C*) or dNy(C) or dNy(CO) 
over (dNy(C^)+dNy(C)+dNy{CO)) averaged over the map. 
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Fig. 9. Channel map of the fraction of column density against the total carbon column density (dNv(C*)/(dNr(C*) + dNy(C)+dNv(CO))) using 
the column density derived from 7’ex= 67 K and assuming rj = 0.5. The velocity bin is 1 kms^*, and the central velocity is shown in each panel. 
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Fig. 10. The fraction of total column density of C^, i.e., N{C^)l{N{C^)+N{C)+N{CO)), for different r}. 


a carbon abundance of 7.9 x 10“^ (iGarnetll 19991) for selected po¬ 
sitions and rj (Table |2]). For a hydrogen density of 10^-10^ cm“^ 
(see above), the depth of the cloud is 0.3-3 pc for 10^^ cm“^. 
The typical projected size of clumps in the observed region is 
30", which corresponds to 7 pc at a distance of 50 kpc. Assum¬ 
ing that the projected size is on the same order as the depth of 
clouds, a smaller 77 would be preferred. The total mass, obtained 
by integrating over the whole observed area, is 2.2 x IO^Mq, 


4.4 X 10^ Mq, and 8.5 x IQ^Mq for 77 = 0.5, 0.2, and 0.1, respec¬ 
tively. 

With dNt, of each species, we estimate the fraction of each 
species in each velocity bin. Figure |9] shows the channel map of 
the column density fraction against the total column density 
of the three species, i.e., dNviC^)/(dNviC^)+dNviC)+dNy(CO)), 
in the case of 77 = 0.5 and the column density derived from 
the constant T^x of 67 K. There are two prominent features here. 
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One is given by dips in the column density fraction at the 
N159 W core (~ 240 hms^'), N159 E (~ 233 kms-‘), and 
southwest of N159 E (~ 235 kms^'). The fraction drops 
to 30%, 50%, and 40%, respectively, and the rest is shared by 
[Ci] (30%, 20%, 20%) and CO (40%, 30%, 40%). On the other 
hand, is dominant at all velocity bins in the region between 
N159 W and E. The second feature is that is clearly domi¬ 
nant at almost all positions in the map for velocities far from the 
line center (< 230 kms“^ and > 240 km s '). The total column 
density depends on the assumed rj, but the fraction of the column 
density of each species is less sensitive to 77 , and this spatial and 
spectral trend is largely independent of the value assumed for rj. 
Eigure [To] shows the fraction of the total against total carbon 
column density, namely N{C^)l{N{C^)+N{C)+N{CO)) for dif¬ 
ferent values of 77 , confirming the spatial trend, and Eig.ITTIshows 
the spectra of the column density fraction averaged over the map, 
which clearly demonstrates that is the dominant species at 
velocities of below 230 km s“' and above 240 km s“'. 


3.4. Estimate of the [Cii] emission coming from ionized gas 


We examine the ionized gas contribution to the [C ii] emission 
using the [N ii] 205 gm emission. Since the [N ii] emission is very 
weak, we convolve the [N ii] emission to a spatial resolution of 
50" and a spectral resolution of 2 kms *. The integrated inten¬ 
sity map (Eig. [T2l l shows an anti-correlation between the [C ii] 
and [N ii] emission, which indicates that most of the emission 
does not come from the ionized gas, which also emits in [N ii]. To 
estimate the contribution of the ionized gas to the [C ii] emission 
quantitatively, we compare the observed [N ii]/[C ii] ratio with 
predictions based on a simple isothermal three-level system, for 
electron temperatures in the range 7000-12000 K. In Table|3| the 
observed ratios (in 50" beams) are listed for the four positions 
indicated in Eig. [121 as well as the average over the blue box area 
in Eig.[T2|and the average over the whole observed area. Eor the 
calculations we used the ele mental abundan ces of C and N for 
EMC H II regions as given by Game td (1 19991). The N abundance 
is consistent with IVermeii & van der Hulsd (l2002h . We assume 
optically thin emission and LTE conditions (high density limit). 
These assumptions provide an upper limit on the contribution of 
the ionized gas to the [Cii] emission for the following reasons. 
When we calculate [Nii] 205/7m/[Cii] 15877 m as a function of 
the electron density, it has a maximum around 50-100 cm“^ and 
reaches its minimum (0.074) in the high density limit. The [N ii] 
20577m/[Cii] 15877 m ratio is much less sensitive to the electron 
density than [Nii] 12277 m/[Cii] 15877 m; the former varies by 
only ~ 65% over the whole electron density range, whereas the 
latter varies by one order of magnitude. Still, if the electron den¬ 
sity is at a range where the predicted [N ii] 205 gm/[C ii] 158 77m 
ratio is higher than 0.074, the [C ii] intensity attributed to the 
ionized gas becomes lower. If the observed [C ii] emission is op¬ 
tically thick (or T ~ 1, see previous section), the optical-depth- 
corrected [C ii] emission would be stronger, so the fraction of 
the ionized gas contribution would decrease. The maximum pos¬ 
sible contribution from the ionized gas to the [C ii] emission is 
40-60% at positions A, B, and D, where the [N ii] emission has 
a peak, with only 20% at position C. Position C is the [C ii] peak, 
which has only very weak corresponding CO and [C i] emission 
(see p7 in Sect. 13. 11 1. Using the [N ii] emission, we exclude the 
ionized gas as a dominant origin of the [C ii] emission at this 
position. The average contribution over the blue box in Eig. [121 
where the [N ii] is detected (at 50" resolution), is 30%, and over 
the whole observed area it is 15% (Table[2l. 


In Eig. [ 13 ] the velocity profiles of [C ii] and [N ii] at these 
four positions are compared. We fit the [N ii] line with a single 
Gaussian, and the fitted central velocity is indicated as a verti¬ 
cal line. Eor positions A and B, the difference of the velocity 
profile between the [N ii] and [C ii] is within the noise, but for 
positions C and D, [N ii] has a higher central velocity, and the 
[C ii] profile of the higher velocity wing matches the [N ii] pro¬ 
file, indicating that this velocity component might be associated 
with the ionized gas. The ratio of [N ii]/[C ii] in the high-velocity 
wing (> 245 km s”') for position C and D is consistent with their 
origin in the ionized gas. 

The column density in the ionized gas can be directly 
estimated by the column d ensity (A(H^)) derived from the 
radio continuum observations. iMartm-Hernandez et al.l d2005h 
derived the emission measure and the electron density at three 
compact peak positions by the 6 cm continuum observations 
usi ng the ATCA. Th e se are the compact H ii regions identified 
by lindebetouw et al.l (l2004l) and match the [C ii] peaks (blue 
asterisks in the CO(4-3) panel of Eig. [TJ p4, pll, and north¬ 
east of pi). The electron density (tie) at these positions is (2.2- 
3.0)xl0^ cm“^, and together with the emission measure it gives 
A(H+) of (3-9 )x102i cm-2, and thus A(C+) of (3-7)xl0i^ cm^^ 
if we assume that all carbon atoms are singly ionized. Since two 
of them are spatially resolved by ATCA having the source size 
of 4-5" and the other one is smaller than the beam size (< 2"), 
the dilution factor in our 20" map is 16-103, which makes the 
estimate of A(C^) = 4 x 10*® cm“^ for pll and northeast of 
pi, and 3 x 10*® cm“^ for p4. They are 2% or less compared to 
the A(C^) derived from the [Cii] emission (Table [2]). However, 
this should be interpreted as a contribution from the dense ion¬ 
ized gas traced preferably by the radio continuum (its intensity is 
proportional to ti^), but not from the diffuse ionized gas. The Ha 
map by the Magellanic Cloud Emission-Line Survey (MCELS) 
shows a different morphology including more ex tended structure 
and additional peak positions (IChen et al.ll2010l) . Our [N ii] map 
does not have a peak at these compact H ii regions, but corre¬ 
lates with Ha, confirming that they instead trace more diffuse 
ionized gas than the radio contin uum because the critical density 
of the [N ii] 205 77 m is 180 cm“^ (IShibailll992l) . Since the critical 
density of the [Cn] 15 8 77 m with the collision partner of elec¬ 
tron is 40 ([Shibailll992h . it also traces more diffuse ionized gas 
rather than dense ionized gas. Thus, the [N ii] 205 77 m emission 
is useful for estimating the ionized gas contribution to the [C ii] 
emission, tracing the same gas. The reason that the [Cii] emis¬ 
sion has peaks at the positions of compact H ii regions should be 
interpreted as a dominant contribution of the PDR gas that sur¬ 
rounds these dense H ii regions and is illuminated by strong the 
UV radiation. 

Theoretical models of Hii region/PDR complexes of solar 
metallicity gas predict contrib utions from the i onized gas to the 
[C n] emission as high as 40% (lAbel et al.l2005h or less than 10% 
(iKaufman et alJUOl^ for most of the physical parameter range. 
Observations toward Galactic PDRs have i n fact provided com- 
parable results. Towards the Carina nebula. lOberst et ^ (1201 Ih 
derive 37% from the [N ii] 205 77m/[C ii] 158 77 m ratio, while us¬ 
ing the correlation of JC n] 1 58 um with [Nn] 12277 m and [Oi] 
63 77 m [Mizutani et al.l (l2004t) find 20%. iGiannini et ^ (l2()00l) 
show a 10-30% contribution from the ionized gas in NGC 2024 
by CLOUDY modeling. In S171 and cr Sco region, the PDR 
contribution to the [C ii] emission is estimated by the PDR mod¬ 
eling, and the ionized gas contribution shows a spatial distribu¬ 
tion with about 80% at positions closer to th e exciting stars, but 
only 20% or 30% at positions farther away (lOkada et al.ll2006l 
120031) . Similarly 10-30% contribution of the ionized gas is re- 
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Table 3. Maximum possible contribution of the ionized gas to the [C ii] 
emission based on the integrated intensity ratio of [N ii]/[C ii]. 


Position 

[Nii]/[Cii] 

Max. ionized 
gas fraction [%] 

A 

0.031 +0.004 

44 + 5 

B 

0.043 + 0.006 

61+8 

C 

0.014 + 0.002 

19 + 3 

D 

0.028 + 0.005 

39 + 6 

(average) 

blue box in Pig.fT2l 

0.022 + 0.002 

30 + 3 

whole region 

0.011+0.001 

15 + 2 
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3.5. Spatially resolved I[cn]/Ii 
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Fig. 13. [Cii] (purple) and [Nii] (black; scaled by 10) spectra at posi¬ 
tions indicated in Fie. 1121 with 2 kms“* spectral resolution. Blue lines 
show the result of a single Gaussian fit to the [N it] emission, and its 
center is indicated as a vertical line. 
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Fig. 12. Integrated (220 - 250 kms^') intensity map of [Nii] (colors), 
convolved to 50" angular resolution, overlayed with the contours of 
[C It] at 20" resolution. Blue asterisks show the positions where spectra 
have been extracted for Fig. The inserted spectrum shows the [N ii] 
emission averaged in the area marked by the blue box (x-axis is Velocity 
[km s“'] and y-axis is Tinb [K]). 


ported in several star-f orming regions in nearby galaxies with al - 
most solar metallic itv (IContursi et al.ll2013l:lKrainer et al.ll2005h . 
iLebouteiller et al.l (12012 ) present that 95% of the [C ii] emission 
arises in PDRs in LMC-Nl IB, except toward the star cluster for 
which as much as 15% could arise in the ionized gas. In M33 
star-forming region s a typical ionized gas contribution is esti¬ 
mated to be 5-50% (iHigdon et al.ll200^ . In principle in a lower 
metallicity environment, the [C ii] emission from the H ii region 
decreases because the size of the H ii region is independent of 
the metallicity and the metal emission line is proportional to the 
abundance, while the column of layer in a PDR stays con¬ 
stant. Therefore, the ionized gas contr ibution in a lower metal¬ 
licity environment s hould be lower (llsrael & MalonevI 1201 It 
I Kaufman et al.l1200^ . However, because of the uncertainties of 
the results cited above, this dependency has not been demon¬ 
strated by observations. 


o 

o 

o 

CM 


o 

LU 

Q 


46 


-69° 48' 



05" 40"’ 20" 


00" 39"" 40" 

RA (J2000) 


The ratio [Cii]//fir has been widely used as a tracer of the 
photoelectric heating efficiency, under the assumption that [C ii] 


Fig. 14. [C ii]//FiR,th rnap (see text for the definition) overlayed with the 
contours of the [C ii] integrated intensity. 


158 pm is the dominant cooling line. Here we derive the spatial 
distribution of this ratio. We obtained the archival data of MIPS 
24/rm (the legacy survey of SAGE), PACS 100/im and 160/im, 
and SPIRE 250/rm (the open key time program HERITAGE); 
applied the color correction; convolved them to a spatial resolu¬ 
tion of 20"; fitted them with a modified blackbody with a single 
temperature and a spectral index of /? = 1.5; and integrated the 
fitted curve over the total range of the thermal emission. Tra¬ 
ditionally, /fir is defined as the 42 pm-122pm flux and /tir is 
defined as the 3/rm-llOOjum flux (iDale & Heloull2()()^ . Here 
we use the direct integration of the fitted total thermal emis¬ 
sion, and express it as /fir,*- We do not calculate /tir because 
we do not have a spectrum of polycyclic aromatic hydrocarbons 
(PAHs) at each position, and its contribution to total infrared flux 
(TIR) can vary in spatially resolved regions. We do not adopt 
/fir either because there is no reason to cut the integration at 
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122/tm. Figure [14] shows the spatial distribution of /[ciii/fpiR,*- 
The ratio is significantly lower at the N159 W and E cores, be¬ 
ing 0.3-0.4%. At the [Cii] peak between N159 W and E (p7), 
7[cii]/7FiR,th = 0.7-0.8%. These values are in the same range as 
the results in the LMC/Nl 1 region by iLebouteiller et ^ (1201 2h 
when considering that TIR, which includes mid-infrared emis¬ 
sions, is 10-35% higher than the to tal thermal emissi on (esti¬ 
mate among the sample of PDRs in lOkada et al.ll207^ except 
for Ced 201, whi ch is excited by a much later star of B9.5). 
Ilsrael et ^ d 19961) present a value of /[cu]/7fir of 1-2% for the 
N159 E and W cores. Their much larger value is due to the defi¬ 
nition of /fir and the lower spatial resolution (the FIR emission 
has a sharper peak at the cores than [C ii]). 

There are three possible causes of the lower /[cu]//FiR,th at the 
cores. The first possibility is a contribution of the [Oi] 63 yum to 
the cooling. Here we do not directly use the archival PACS spec¬ 
troscopic data because they do not cover our whole observed 
region, but the comparison of the [O i] 63 yum and [C ii] 158 yum 
emissions within their observed area indicates that the intensity 
of [Oi] 63yumis 30-100% of the [Cii] 158//m intensity. We plan 
to observe velocity resolved [O i] 63 yum with GREAT at our next 
southern flight opportunity and carefully compare the line pro¬ 
file. The second possibility is that the photoelectric heating effi¬ 
ciency is actually low at these cores. The photoelectric heating 
efficiency is a function of the charging parameter y = jug, 

where Go is UV radiati on held, T is the gas te mperature, and Ug 
is the electron density (iBakes & Tielensl|l994h . When y is high, 
dust grains are positively charged. Thus, the photoelectric heat¬ 
ing efficiency becomes lower because a higher energy is needed 
to further ionize them. If there is a strong UV source in the core, 
y in the core is high and the heating efficiency could be lower 
than outside. This may be likely because both N159 E and W 
possess compact H ii regions (Eig.[T]CO(4-3) panel). Correlation 
with the fraction of positively ionized PAHs would give more 
detailed examination of this possibility (lOkada et al.ll2013l) . but 
there are only sparse Spitzer/IRS observations in this region. The 
third possibility is that later-type heating sources contribute to 
fFiR.th but not to the photoelectric heating because of the lack of 
photons with the energy of > 6 eV. 

4. Summary 

We mapped a 4' x (3'^') region in the N159 star-forming re¬ 
gion in the EMC in the [Cii] 158//m and [Nii] 205//m lines 
with the GREAT instrument on board SOEIA, as well as CO(3- 
2), (4-3), (6-5), '^CO(3-2), and [Ci] ^Pr^Po and 3P2-^Pi with 
APEX, and study their velocity profiles in detail. The variation 
of the velocity structure over the observed area is large, and the 
[Cii] emission shows a significantly different velocity profile 
from the other lines at most positions. Quantitatively, we esti¬ 
mate the fraction of the [C ii] emission that cannot be attributed 
to the material traced by the CO line profiles to be 20% around 
the CO cores and up to 50% in the area between the cores. We de¬ 
rived the column density of C^, C, and CO in each velocity bin, 
and And that the relative contribution from C^ dominates at the 
velocity range far from the line center, and the area between the 
CO cores. The origin for this [C ii] emission, at velocities sub¬ 
stantially off the velocities of the dense molecular cores traced by 
CO and [Ci], may be due to gas ablating from the dense cloud 
cores or to additional gas not immediately associated with the 
dense cores. 

To estimate the contribution of the ionized gas to the [C ii] 
emission, we use the ratio with the [N ii] 205 //m line and show 
that the ionized gas can contribute up to 19% to the [Cii] emis¬ 


sion at its peak position, less than 30% in the area where [N ii] 
is detected, and less than 15% if averaged over the whole area 
observed. By comparing the peak velocity of the [C ii] and [N ii] 
emissions, we provide evidence that a high-velocity wing in the 
[C ii] emission is likely due to ionized gas. Our study emphasizes 
that it is essential to resolve the velocity profile in order to distin¬ 
guish between different cloud components and correctly derive 
the physical properties of individual components. 
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